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SUMMARY 


A potential-flow panel method has been modified to calculate the effects 
of a rotor wake on the time-averaged surface pressure and velocity distribu- 
tions on a helicopter fuselage. The rotor-induced velocities are calculated 
by using a vortex-tube wake model. The calculated pressure distributions are 
found to compare well with experimental data obtained from tests of a wind- 
tunnel model. 


INTRODUCTION 

The fuselage of a helicopter operates in a flow field which is composed 
of a combination of its forward flight velocity and the downwash induced by the 
rotor system. This combination can range from only the rotor-induced field at 
hover to very little rotor downwash effect at high speeds. Analysis of this 
flow field at any condition from hover to high-speed flight is important in 
designing the fuselage and its various components (such as wings and tails) . 

Potential-flow panel methods (ref. 1) , which use source and/or doublet 
panels to model rather complex vehicles, have been developed to analyze the 
flow field of fixed-wing aircraft. Much work (refs. 2 to 12) has been done in 
the helicopter community to analyze the induced flow of a rotor system. These 
methods range from simple actuator disk theory to unsteady lifting-surface 
methods . 

Even though the dynamic nature of the rotor wake is important in analyzing 
certain problems such as vibration, many design details can be analyzed by 
using a time-averaged downwash field. Therefore a computer code has been 
developed to combine a vortex-tube rotor-wake theory with an incompressible, 
potential-flow panel method. This method calculates both on-body and off-body 
velocities , fuselage surface pressures , and total loads . 

This paper includes a discussion of the basic theories, which are incor- 
porated in the present computer code, with a general description of the com- 
puter code and the experimental test, which was used as a basis for validation. 
Experimental fuselage surface pressure data were obtained from a helicopter 
wind-tunnel model at an advance ratio of 0.05. Comparisons of experimental 
data and the analytical method are presented for thrust coefficients of 0.0034, 
0.0050, 0.0066, and 0.0082. 


SYMBOLS 

Units used for physical quantities defined in this paper are given in the 
International System of Units (SI) . The positive senses of parameters are 
shown in figure 1. 
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coning angle, deg 

first-harmonic rotor longitudinal flapping angle, deg 
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rotor-wake cross-section area, m 
lateral cyclic control, deg 

first-harmonic rotor lateral flapping angle, deg 
longitudinal cyclic control, deg 
,C 8 constants defined in equation (37) 
pressure coefficient 

rotor thrust coefficient, t r/[p ^ r > 2|rR 5 

= Axial distance of point from axis of vortex ring/R 

= Radial distance of point from plane of vortex ring/R 

= Shortest distance from point to vortex ring/R 

= Longest distance from point to vortex ring/R 

elliptic integral of second kind 

height of hub above gimbal pivot point, cm 

= Height of fuselage/R 

rotor axial force, N 

elliptic integral of first kind 

tip loss factor, Effective rotor radius/R 

nonrial vector 

ellipse power or slipstream contraction factor 

static pressure, Pa 

total pressure. Pa 

point designation (figs. 2 and 3) 

dynamic pressure, O.SpV^ , Pa, or point designation (fig. 2) 
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= Radial polar coordinate/R or Wake radius/R or control point 
distribution (see fig. 2) 

rotor radius or vortex-ring radius, m 

= Root cutout/R 

radial distance of point from axis of vortex ring, m 

incremental surface 

rotor thrust, N 

velocity vector, m/sec 

momentum induced velocity, m/sec 

rotor tip speed, flR, m/sec 

component of V, m/sec 

= Width of fuselage/R 

Cartesian coordinates in body axis 

coordinates in TPP system 

distance from moment reference center of gimbal pivot point, cm 
rotor side force 
— Body camber/R 

angle of attack, deg (see fig. 1(a)) 
angle of sideslip, deg 

rotor-shaft tilt angle, deg (see fig. 1(b)) 
vortex strength 

rotor inflow ratio, sin a - ^xjj V T 

rotor advance ratio, cos a/V T 

air density, kg/m*^ 
source strength 

= ( d 2 " d l)/( d 2 + d l) 

azimuthal polar coordinate, deg, or velocity potential 


3 



■ 


$ total -velocity potential 

X wake skew angle measured from vertical axis, deg 

\p blade azimuth position, deg, or wake-area contraction ratio, A/A q 

(see fig. 1(b)) ; 

¥ stream function 

Q rotor rotational speed, rad/sec (see fig. 1(b)) 

Subscripts : 
f fuselage 

i local condition 

n normal 

o at rotor disk 

p perturbation due to fuselage 

R rotor 

s surface 

t tangential 

w wake 

00 free stream or fully contracted wake 

Notation : 

eff effective 

HP hub plane 

NFP no feathering plane 

rpm revolutions per minute 

TPP tip path plane 

j | vector magnitude 

A difference due to rotor wake 
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gradient 


( ) vector 


DISCUSSION OF THEORY 

The numerical analysis method presented in this paper is based upon the 
panel method of reference 1 and an extension of the vortex-tube rotor-wake 
model derived from references 3 and 4. 


Panel Method Theory 

The panel method is derived from the conditions that the velocity poten- 
tial $ must satisfy Laplace’s equation 

V 2 $ = 0 (1) 

as well as two boundary conditions. The first boundary condition is the 
requirement of tangential flow at the impermeable body surface and can be 
written as 



The second boundary condition is that $ must approach the nondisturbed free- 
stream potential at infinity. The velocity potential $ is comprised of the 
free-stream potential a nonuniform onset or rotor-wake potential 4> w , and 

a disturbance potential due to the fuselage Therefore <$> can be 

expressed as 

4 = ^oo + <f> w + <£p (3 > 

The gradient of the uniform onset potential is given in reference 1 as 

V4>oo - -V^ (4) 

and the gradient of the nonuniform onset potential is defined such that 

V<J> = -Av (P) (5) 

w 

where Av is dependent upon the position of point P. The disturbance poten- 
tial is given as a source density distribution over the body surface such that 


V p) 



(q) 

(P,q) 


( 6 ) 
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where P is a point in space and q is a point on the surface which is a dis- 
tance r from point P. (See fig. 2.) The function a (q) is the source 
density distribution. 

Since equations (4) to (6) satisfy Laplace's equation and the infinity 
boundary condition, the source function G can be calculated by using the 
boundary condition of tangential flow at the body surface (eq. (2)). For 
point P on the body surface, the normal derivative of equation (6) is shown 
in reference 13 to be 
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The normal derivatives of <t> and 6 are 

T oo w 


and 


9<t> 0 

3n 
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(8) 


( 9 ) 


respectively. Thus, combining equations (7), (8), and (9) gives 


2TT G (P) 



3_ 

3n 


r (P , q) 


G (q) ds = -n 



( 10 ) 


where the right side of equation (10) is the normal onset velocity (uniform 
plus nonuniform) and the left side of the equation is the normal perturbation 
velocity. 

In panel methods equation (10) is typically solved by approximating a 
body surface with discrete quadrilateral panels. Thus, equation (10) can be 
expressed as a matrix equation of the form 



where N is the total number of panels and . is a matrix of aerodynamic 
influence coefficients which is a function of geometry only. (See ref. 13 for 
details.) Therefore equation (10) has been replaced by a set of linear 
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algebraic equations. This system of equations is solved by standard numerical 
matrix techniques to obtain the source strengths CJj . 

The velocity vector at each panel control point is obtained by using the 
source strengths and aerodynamic influence coefficients. These final veloci- 
ties are given by 

N 

v± = Yj ^ij CT j + Vo + Av± (i = 1, 2, . . N) (12) 

j=l 

The panel pressure coefficient c p,i is defined as 

C_ • = ~ ?£°. (13) 

Pfl ^co 


Since the static pressure p^ at the panel center can be expressed as the 
difference between the total and the dynamic pressure at panel i 




„ p t,i ~ qi - p 0 

<3co 


(14) 


Also, the total pressure p t ^ can be expressed as the sum of the uniform onset 
flow total pressure p t plus a local total-pressure variation Ap t ^ due to 
the rotor wake. Therefore 


r . = p t,co " ^i - P co + Ap t,i 

P,x ' ' " 


(15a) 


or since p 

t f oo 


q + p f 

CO oo 


~P/1 


= %o + Poo - qi - Poo + A Pt,i 


The free-stream static pressure p^ is eliminated , thus giving 


(15b) 


“Pii 


= - qi + A Pt,i 

‘Sco 


(15c) 
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gives 


and by replacing with 


o.spv ^ 2 


(0.5)plvJ 2 - (0. 5) p | V ± | 2 + Ap t#i 
(0.5)p|vJ 2 < 3oo. 


which reduces to 


C 


P/i 


1 



%> 


(15d) 


(15e) 


The term Ap t j_ will be discussed in the section concerned with the rotor- 
wake modeling! Body forces and moments are calculated by summing the pressures 
over the surface of the body. 


Vortex-Tube Theory 

The effect of the rotor flow field is added to the panel method solution 
of the body aerodynamics as an onset flow disturbance Av^ (eqs. (11) and (12)) 
and as a total-pressure variation Ap tf i/q^ (eq. (15e) ) at each panel control 
point. The onset flow disturbance is determined from the vortex-tube method, 
and the total -pres sure variation is determined from momentum theory. 

The vortex-tube method of wake representation which was used is a varia- 
tion of the method described in references 3 and 4. In the vortex-tube method, 
the wake is defined by a series of concentric vortex tubes each of which con- 
sists of a finite but large number of vortex rings. (See fig. 3.) 

The stream function at a point relative to a vortex ring is shown in ref- 
erence 14 to be 


Y(P) = -?l( d l + d 2 ) [ K(T) ‘ 


(16) 


where T is the strength of the vortex ring, R is its radius, d-^ and d 2 
are the minimum and maximum distances, respectively, from point P to the 
ring (see fig. 3), and K(t) and E(T) are complete elliptic integrals of the 
first and second kinds, respectively. The parameter T is defined as 


d 2 ~ d l 
d 2 + d 1 


(17) 
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The axial component of velocity of the stream function 'V is 
1 


v n = - 


Rp dRp 


where Rp is the radial distance of point P from the vortex-ring axis, 
reference 6, equation (18) is expressed as 


( 18 ) 


In 


V_ = _£ (AB + CDF) 

n 2TTR-D 


(19) 


where 


A = K(T) - E(T) 


B = 


d~ - 1 d- + 1 




C - d]_ + d 2 


D = TE(T) 

1 - T 2 


( 20 ) 

(21) 

( 22 ) 

(23) 


F = 1 - 


‘1-D 
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- d x d 2 _ (1 + d r )d x 2 - (1 - d r )d 2 2 

(24) 

ir 2 

2d r d l d 2 


1/2 


h J 


(25) 

ol 

1/2 


1 ) 2 J 


(26) 


where d a is the nondimensional axial distance from the plane of the ring to 
point P and d r is the nondimensional radial distance from the axis of the 
ring to point P. 


The axial component of velocity for a vortex tube can be expressed as 


V » - 2 (Affl ♦ CjDjFj) 

j=l 


(27) 
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where K is the total number of vortex rings. The vortex-ring circulation T 
is obtained by requiring that V n be equal to the average induced velocity at 
the center of the rotor disk with uniform loading. 

Reference 4 shows that nonuniform disk loading can be represented by a 
series of concentric vortex tubes. Loading is either added (positive vorticity) 
or subtracted (negative vorticity) from the tip vortex-tube strength to vary 
the loading across the rotor disk. 

In the present method , wake contraction was added by varying the vortex- 
ring radii according to empirical relationships given in reference 15. The 
ratio of the radius of the final contracted wake to the rotor radius is 
expressed as 

r ( -Cm* 58 . 77\ 

— = 0.707 + 0.1418^1.0 - e J (28) 

r o 

or in terms of area 



The area ratio at a given axial location z R is expressed as 



z p *N 

\\) + (1.0 - \p)e 


(30) 


where N is a slipstream contraction factor based on experimental data given 
in reference 15. 

Swirl is introduced as a tangential component of induced velocity since 
the lift vector of a blade element is in reality not perpendicular to the disk 
plane but is perpendicular to the local resultant velocity (ref. 15) . There- 
fore, the local tangential velocity V t can be calculated by the equation 



where Vj is the average induced downwash at the rotor disk calculated from 
momentum theory. 
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The pressure coefficient used in the panel method was given as equa- 
tion (15e) which includes a total -pressure term. The derivation of the total- 
pressure term is given in the appendix and is expressed as 


A p t ,i 

qoo 


v + Av 

oo 


2 




(32) 


Equation (32) is evaluated at each panel control point. 


COMPUTER CODE DESCRIPTION 

Analytical results were obtained by using the panel method of reference 1 
coupled with the vortex-tube method as previously discussed. In the computer 
code architecture, the wake calculation program was added as a subprogram to 
the panel method program. The wake program is called after all the body 
geometry has been input and the control points (panel centroids) have been 
calculated. The vortex program then calculates rotor-induced velocities and 
total-pressure variations at the panel control points and returns these values 
to the main program. Source distributions are calculated based upon the total 
onset flow, and velocity distributions are calculated based upon these source 
distributions and the onset flow. The pressure distributions are then calcu- 
lated based upon the velocity distribution and the differential total pressure 
as calculated in the vortex-tube program. 


Modifications to the Panel Method Program 

Modifications to the panel method program included the addition of arrays 
to store the nonuniform onset flows and total -pressure values, addition of the 
nonuniform onset flow term to equations (11) and (12) , and the addition of the 
total-pressure variation term to equation (15e) . 


Vortex-Tube Program 

The vortex-tube program requires as input: (1) Coordinates of panel cen- 

ters; (2) rotor geometry including hub position, tip-path-plane angle of attack, 
fuselage angle of attack, disk radius, and root cutout; and (3) performance 
parameters including C T , V T , tip loss factor, and radial blade loading. 

Initial calculations determine general rotor-wake parameters which define 
the wake geometry. These parameters include the average induced velocity, the 
contraction ratio, the wake skew angle, and the static-pressure differential at 
the rotor disk. 
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The average induced velocity is calculated by using the momentum theory. 
For hover 


V I - V T 


(°. 


5 * C T A eff7 


U/2 


(33) 


where A e ff = l.OyMKp 2 - R c 2 ) is a correction for the effective rotor area, 
K t is a tip loss factor, and R c is the root cutout. 


For forward flight 


V z = 0.5 


Vip * Cip A eff/ 


^(y 2 + A 2 ) 


1/2 


(34) 


which is solved in an iterative fashion since X is a function of Vj. Equa- 
tion (33) is used as an initial estimate of Vj to calculate the inflow 
ratio X. The inflow ratio and induced velocity are alternately calculated 
until Vj converges to a solution. 

Initial ob servations of the skew angle measured from the wind-tunnel test 
(which is discussed later) indicated that the skew angle is well approximated 
by the expression 


X = tan 



(35) 


where ip is the contraction ratio and is calculated by using equation (29) . 


The following calculations are repeated for each panel control point. 
After the coordinates of the control point are rotated and translated to the 
rotor tip-path-plane (TPP) coordinate system, the axial induced velocity is 
calculated by using 1 to 15 concentric vortex tubes (depending upon the radial 
variation of disk loading) . Each tube consists of 101 discrete vortex rings, 
the radii of which vary with the expression 


fi = [iJj + (1.0 - ip)e ZRN ] 


(36) 


where z R is the axial distance of the ring beneath the rotor and N is the 
empirical slipstream contraction factor. Therefore, the total axial induced 
velocity at one control point is the sum of the contributions of 101 vortex 
rings times the number of vortex tubes. 

After the induced axial velocity is obtained, the tangential velocity 
(swirl) is calculated by using equation (31) . The axial and tangential veloci- 
ties are then rotated into the fuselage Cartesian axis system. The local 
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total -pressure variation is calculated by using equation (32) . These rotor- 
induced velocities and total-pressure values are then returned to the panel 
method program. These rotor computations add approximately 5 percent to the 
running time of the panel method program. 


EXPERIMENTAL DATA 

Validation of new analytical methods generally requires experimental data 
as a basis for comparison. Although reference 16 presents surface pressures 
of three helicopter models with a 1-meter-diameter , two-bladed rotor, the 
pressure orifice locations were considered too sparse to allow for accurate 
comparisons with the analytical method. Therefore a wind-tunnel investigation 
using a helicopter model with more closely spaced pressure orifices was con- 
ducted in the Langley V/STOL Tunnel. A photograph of the model in the test 
section is presented as figure 4. 


Model Description 

The rotor system used in this investigation consisted of a 3.15-meter- 
diameter, four-bladed rotor. The hub was fully articulated with flapping and 
lagging hinges coincident at 4.8-percent radius. The rotor blades had 
-8° twist, an untapered planform, and square tips. Details of the rotor system 
may be found in table I. 

The fuselage shape for this test is mathematically defined. At a given 
fuselage station x, the cross-section y- and z~coordinates are defined by the 
local fuselage height H, width W, camber line ZO, and elliptical power N 
as follows. 


The parameters H, W, ZO, and N are obtained by applying the function 


f (x) 


C 6 



+ C. 


'x + C 3 \ C 5 
. C 4 / 


i/c 8 


(37) 


with separate sets of constants to Cg for each of the four cross-section 

parameters H, W, ZO, and N. As shown in figure 5, the fuselage is divided 
into four regions and the pylon into two regions with a set of constants for 
each region. These constants are given in table II. 

By using the three parameters H, W, and N as obtained from equa- 
tion (37) , the cross section can be defined by the polar coordinates r and (J) 
such that for a given $ 


r(cf>) = 


(0.5H + 0.5W) 


N 


(0.5H sin) N + (0.5W cos) 


N 


1/2 


(38) 
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( 39 ) 


The body Cartesian coordinates are then obtained from the relationships 
y = r(<f)) sin <f) 

z = r (<#)) cos cj) + ZO (40) 


Test Conditions 

The experimental data presented in this paper were obtained for an advance 
ratio of 0.05 at the test conditions as given in table III. The wind-tunnel 
test section , which measures 4.42 meters by 6.63 meters, was configured with 
the walls and ceiling removed. Angle of attack and test-section dynamic pres- 
sure were corrected for wall effects by using methods described in reference 17. 
All tests were conducted with the moment reference center of the model on the 
center line of the test section. This position corresponds to a rotor height 
of 1.84 times the rotor radius. The pressure data were averaged for 50 samples 
over 5 seconds of sampling time. The data used in this paper were obtained from 
reference 18. 


PRESENTATION OF RESULTS 

The results obtained from the computer code and the experimental tests are 
presented in the following figures: 


Figure 


Isolated rotor-induced velocities in y = 0 plane for y = 0.05 ... 6 

Calculated x- and z-components of total velocity on left side of 

fuselage for y = 0.05 7 

Pressure distribution along top center line of fuselage 

for y = 0.05 8 

Lateral pressure distributions at four fuselage stations 

for y = 0.05 9 

Separation and swirl effects on fuselage cross-section pressure 

distribution 10 


DISCUSSION OF RESULTS 

Results obtained from the analytical method are presented in the form of 
rotor-wake onset velocity distributions, total-surface velocity distributions, 
and surface pressure distributions. These three outputs are useful in analyzing 
the aerodynamics of a helicopter in different ways. 

Figure 6 presents the induced velocity vectors of the isolated rotor in a 
vertical plane along the center line of the hub (i.e., y R = 0) at four thrust 
levels. Evident in this figure is the contraction and skew angle of the rotor 
wake as calculated by the rotor-wake segment of the program. As would be 
expected, the skew angle decreases with increasing thrust as the momentum 
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inflow velocity increases. The skew angle ranges from approximately 40° at the 
low thrust to 25° at the high thrust. Full wake contraction occurs at approxi- 
mately 50-percent radius below the rotor disk. 

The inflow variation at the rotor disk is also evident with zero inflow at 
the center and maximum downwash at approximately 80-percent radius. There is 
also a substantial difference in inflow levels for the front and aft portions 
of the rotor. These variations decreased downstream of the rotor. 

In figure 7 the x- and z-components of the total-velocity vectors on the 
left side of the fuselage are plotted. These results are obtained from the 
panel method as rotor-wake effects on the fuselage. Wake boundaries are dis- 
cernible on the nose and tail cone as large gradients in the vertical component 
of the total velocity. These total velocities follow the same general pattern 
of the isolated rotor-induced velocities. Longitudinal variations of vertical 
velocity are particularly evident on the tail cone. The effect of the root 
cutout is seen as a stagnation region below and slightly behind the hub loca- 
tion at 35-percent-body length. 

Figures 8 and 9 present comparisons of surface -pres sure coefficients for 
experimental data and results obtained from the analytical method. Although 
the pressure distribution on the entire fuselage surface is not shown, these 
figures illustrate the character of the longitudinal and lateral pressure 
distributions . 

In figure 8 the pressure coefficients and geometry of the top center line 
of the fuselage are given. Since the top center line would be considered the 
stagnation line for a rotor downwash, the large positive values of C p should 
be proportional to the square of the downwash velocity. Very evident in this 
figure is the nonuniformity of the longitudinal downwash distribution as is 
seen in the velocity vectors of figures 6 and 7. 

Even though the overall character of the experimental and analytical dis- 
tributions match, several discrepancies are evident. For the low- thrust case 
(fig. 8(a)), the wake impingement location (point of maximum C p ) is not 
coincident (x/R ^ 0.25) because the analytical model assumes a constant skew 
angle and the actual wake is curved. 

Other discrepancies in figure 8(a) include the overall levels of the down- 
wash with the analytical method underpredicting the pressure coefficients and 
the seemingly higher experimental downwash from the front part of the rotor 
disk as compared with the aft portion. These discrepancies are not evident at 
the three higher thrust levels, although the analytical method does slightly 
underpredict the pressures on the aft portion of the fuselage. This underpre- 
diction may at least be partly due to the probable separation of the rotor 
pylon which is not modeled in the analytical method. 

Experimentally, the forward location of impingement varied from x/R = 0.25 
for C T = 0.0033 to x/R = 0.15 for Crp = 0.0082, which corresponds to skew 
angles of approximately 40° to 25°, respectively. The aft impingement point 
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cannot be discerned from the experimental data as it occurs slightly aft of the 
last pressure orifices. 

Figure 9 gives the pressure distribution, experimental orifice location, 
and panel centroid locations at four fuselage stations and four thrust levels. 
These four stations are x/R = 0.20, 0.30, 1.34, and 1.53. Two are ahead of 
the rotor hub and two are behind. The experimental data, especially at the 
higher values of C^, in these figures generally follow an M-pattern with large 
negative values of Cp where the downwash accelerates around the cross-section 
corners and large positive values of Cp at the top (z/R maximum) and bottom 
(z/R minimum) stagnation points. The analytical data have this same 
M-character but quantitatively differ at two regions. The analytical method 
seems to underpredict the flow acceleration around the top corner (z/R posi- 
tive) and to overpredict the stagnation pressure on the bottom side. The 
reason for the first discrepancy may be the coarseness of the paneling at the 
corners and the fact that the wake program generally underpredicts the downwash 
levels. The second discrepancy is attributed to flow separating at the lower 
corner of the cross section as illustrated in figure 10. As mentioned pre- 
viously, a separation model was not included in this panel method program. 

Also evident in figures 9 and 10 is the effect of swirl. A crossflow 
on a fuselage section would decrease the flow acceleration at the corners on 
one side and increase the acceleration on the opposite side; for example, 

C T = 0.0082. At x/R = 0.20 (see fig. 9 (m) ) , the difference in the experi- 
mental pressures at the top corner for the right and left side is 4.75 with 
the left side being higher (more negative Cp) . The analytical results again 
underpredict these corner pressures. Behind the rotor at x/R = 1.53 
(fig. 9(p)), the pressure orifice spacing was not fine enough to indicate the 
quantitative difference of the right and left side, but the right-side pressure 
distribution indicates a higher acceleration. Therefore the crossflow right 
to left forward of the hub and left to right aft of the hub indicates a defi- 
nite swirl in the direction of rotor rotation as modeled by the wake program. 


CONCLUSIONS 

The surface pressure and velocity distributions about a helicopter fuse- 
lage immersed in a rotor wake have been calculated by using a potential-flow 
panel method having a vortex-tube rotor model. The pressure distributions have 
been compared with experimental results obtained from a 3.14-meter-diameter 
rotor wind-tunnel model with a general helicopter fuselage. The conclusions 
are summarized as follows : 

1. The measured and predicted pressure on the fuselage at low rotor-wake 
skew angles varies significantly both longitudinally and laterally as a result 
of the influence of the rotor wake. 

2. The panel-method — vortex-tube combination, which accounts for skew 
angle, contraction, and total -pres sure variations, can adequately calculate the 
character of the time-averaged pressure and velocity distributions if the tube 
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geometry (skew angle and contraction) is modeled correctly , the total-pressure 
corrections are accounted for, and the swirl velocities are added. 

3. Increased accuracy would be obtained with the inclusion of a flow- 
separation technique and a curved rotor-wake model. 


Langley Research Center 

National Aeronautics and Space Administration 
Hampton, VA 23665 
April 17, 1980 
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APPENDIX 


CALCULATION OF TOTAL PRESSURE IN THE ROTOR WAKE 

Momentum theory in hover considers the rotor as a disk as shown in the 
following sketch: 




Location 0 is at infinity; location 1 is immediately before the disk; loca- 
tion 2 is immediately following the disk; and location 3 is at the point of 
full contraction. The average induced velocity through the disk is Vj , and 
the fully contracted wake velocity is w. The static pressures are p Q , 

P 2 , and p^ at locations 0, 1, 2, and 3, respectively. 


Bernoulli's equation for steady , inviscid, incompressible flow with negli- 
gible gravitational force and no work done is 


p t = p + 0.5pV 2 


(Al) 


where p t is the total pressure, p 
sity, and V is the local velocity, 
tions and assuming that p t Q - Pt,l 

p 0 = Pi + 0.5PV-J- 2 


is the static pressure, p is the den- 
Applying equation (Al) to the four loca- 
and p t f 2 = Pt,3 gives 


(A2) 


and 

P 3 + 0.5pw 2 = p 2 + 0.5PV.J. 2 


(A3) 


The thrust on the disk can be expressed as 


T = ( P 2 " P l) A 
where A is the disk area. 


(A4) 


Subtracting equation (A2) from equation (A3) and solving for f p 2 
gives v 

p 2 - Pi = p 3 - p 0 + 0 . 5pw 2 



(A5) 
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APPENDIX 


However P 3 “ Po = if the assumption is made that the static pressure in 
the fully contracted wake is equal to the ambient static pressure. By making 
this assumption and combining equations (A4) and (A5) , the rotor thrust is 


expressed as 

T = 0.5pw 2 A (A 6 ) 

From momentum theory 

T = pAVjW (A7) 

Equating equations (A 6 ) and (A7) and solving for w gives 

w = 2.0V! (A 8 ) 

Using equation (A 8 ) in equation (A5) gives 

P 2 - Pl = 2.0PV ;,. 2 (A9) 

From equation (A3) , assuming that P 3 = p 0 

P 2 - Pp = 0.5pw 2 - 0.5 PVl 2 (A10) 

Again , by substituting w = 2.0Vj, equation (A10) becomes 


P2 " Pp " 1- 5 P V ! 2 


(All) 


or 


P2 " Po = 0-75(p 2 - p x ) (A12) 

Therefore the static-pressure rise immediately behind the disk is three-fourths 
of the static-pressure rise across the disk. By applying the definition 


C,p (A13) 

pAV T 2 

to equation (A4) , p 2 - Pq becomes 

P 2 - Pi = C T pV T 2 (A14) 

Substituting equation (A14) into equation (A12) gives 

P 2 “ P 0 = 0.75C T pV T 2 (A15) 
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APPENDIX 


or by nondimensionalizing by O.SpV^ and making the approximation y = Voo/Vt, 
equation (A15) becomes 


P2 ~ PQ _ 1-5C t 

o • Spv ^ 2 y 2 


(A16) 


which expresses the difference in static pressure across the wake boundary. 

The difference in total pressure inside the rotor wake at a point i and 
free stream can be expressed as 


A P t ,l _ A Pi + A qi 

^ <Joo ^00 


(A17) 


By assuming that the static-pressure variation at point 2 in forward flight is 
the same as hover and that the static-pressure variation goes to zero for the 
fully contracted wake, p can be expressed as 





(A18) 


where A^ is the wake cross-section area at point i. 

The difference in dynamic pressure can be expressed as 


Aqj = [Vqo + AVj 
^ V 2 

v oo 


(A19) 


where Av^ is the resultant of the axial and tangential wake velocities. 
Therefore 


APt,i 

3co 


V + AV, | 2 

oo 1- 

I v „| 2 



(A20) 
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TABLE I.- MODEL GEOMETRY 


Fuselage : 

Moment reference center: 
x 

Y 

z 

Length 

Rotor: 

Hub coordinates : 

x 

Y 

z 

Number of blades . . . 

Y f deg 


Root cutout 

Chord, m 

Radius , R, m .... 
Twist, deg 

2 

Flapping inertia, kg-m 

Solidity 

Airfoil section . . . 


0.690R 

OR 

OR 

2. OR 

0.690R 

OR 

0.274R 

4 

2 

0.009R 

0.034R 

0.090R 

0.20R 

0.108 

1.574 

- 8.0 

0.653 

0.0871 

NASA RC-10-(B) M002 
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View from left 



View from rear 


(b) Rotor system. 


Figure 1.- Concluded. 








^ L-78-4383 

Figure 4.- Wind-tunnel model installed in test section of Langley V/STOL Tunnel. 
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(c) C™ = 0 . 0066 . 


Figure 6.- Continued 
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(a) C T = 0.0034. 



Figure 7.- Calculated x- and z-components of total velocity on left side of fuselage for y = 0.05. 
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Figure 8.- Pressure distribution along 


0 . 0033 . 


center line of fuselage for y 
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Figure 8.- Concluded 
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Figure 9.- Continued. 
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Figure 10.- Separation and swirl effects on fuselage cross-section 

pressure distribution. 
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